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Abstract 

We analyze the crossover from the Bardeen-Cooper-Schrieffer (BCS) state of weakly bound 
Fermi pairs to the Bose-Einstein condensate (BEC) of molecular dimers for a Fermi gas made of 
neutral atoms in three hyperfine states with a SU(3) invariant attractive interaction. By solving 
the extended BCS equations for the total number of particles and the pairing gap, we calculate at 
zero temperature the pairing gap, the population imbalance, the condensate fraction and the first 
sound velocity of the uniform system as a function of the interaction strength in both three and 
two dimensions. Contrary to the three-dimensional case, in two dimensions the condensate fraction 
approaches the value 1 only for an extremely large interaction strength and, moreover, the sound 
velocity gives a clear signature of the disappearance of one of the three hyperfine components. 
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I. INTRODUCTION 



In the last years degenerate ultracold gases made of bosonic or fermionic atoms have 
been the subject of intense experimental and theoretical research [l, 2|. Among the several 
hot topics recently investigated, let us remind the expansion of a Fermi superffuid in the 
crossover from the Bardeen-Cooper-Schrieffer (BCS) state of weakly bound Fermi pairs to 
the Bose-Einstein condensate (BEC) of molecular dimers 

aa 

; the surface effects in the 

unitary Fermi gas j^]; the localization of matter waves in optical lattices the transition 
to quantum turbolence in finite-size superfluids [7|, |8|. 

Very recently degenerate three-component gases have been experimentally realized using 
the three lowest hyperfine states of 6 Li js), [lo|. At high magnetic fields the scattering lengths 
of this three-component system are very close each other and the system is approximately 
SU(3) invariant. Moreover, it has been theoretically predicted that good SU(N) invariance 
with N < 10) can be reached with ultracold alkaline-earth atoms (e.g. with 87 Sr atoms) 
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13). In the past various authors 
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171 ] have considered the BCS regime of a fermionic 



gas with SU(3) symmetry. In the last years He, Jin and Zhang [18 



19[ have investigated the full BCS-BEC crossover 2CH24) of this system at zero and finite 



and Ozawa and Baym 



temperature in three-dimensional space. Recently we have calculated the condensate fraction 
and the population imbalance for this three-component quantum gas both in the three- 



dimensional case and in the two-dimensional one 



18 



19 



25) . In this paper we review the extended 



3EC theory 2lM23) for an atomic gas with three-component fermions at zero temperature 



25) but without invoking functional integration. We obtain the chemical potential, 



the energy gap, the number densities and the condensate fraction as a function of the 
adimensional interaction strength. Finally, we calculate also the first sound velocity of the 
system both in three and two dimensions. In two dimensions we find that the sound velocity 
shows a kink at the critical strength (scaled binding energy) where one of the three hyperfine 
components goes to zero. 

The Lagrangian density of a dilute and ultracold three-component uniform Fermi gas of 
neutral atoms is given by 

ih di + 2m ^ + ^ J ^ a ~ 9 [^R^G^G^R + ^r^b^b^r + ^g^b^b^g) , 

(1) 

where ip a (r,t) is the field operator that destroys a fermion of component a in the position 



r at time t, while ip£{ r ) creates a fermion of component a in r at time t. To mimic QCD 
the three components are thought as three colors: red (R), green (G) and blue (B). The 
attractive inter-atomic interaction is described by a contact pseudo-potential of strength g 
(g < 0). The average total number of fermions is given by 

N = E (2) 

a=R,G,B J 

where (• • •) is the ground-state average. Note that N is fixed by the chemical potential a 

nn 

which appears in Eq. (JT]). As stressed in Refs. [18|, [19J, by fixing only the total chemical 
potential fi (or equivalently only the total number of atoms N) the Lagrangian ([1]) is invariant 
under global SU(3) rotations of the species. 



II. EXTENDED BCS EQUATIONS 



At zero temperature the attractive interaction leads to pairing of fermions which breaks 



the SU(3) symmetry but only two colors are paired and one is left unpaired llSl, l25|] • We 



assume, without loss of generality 
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25], that the red and green particles are paired 



and the blue are not paired. The interacting terms can be then treated within the minimal 
mean-field BCS approximation, i.e. neglecting the Hartree terms while the pairing gap 



A = g(4> G {T,t)$ R (r,t)) 



(3) 



between red and green fermions is the key quantity. In this way the mean-field Lagrangian 
density becomes 



C 



mf 



( d h 2 \ 

- r n r \ / 



(4) 



a=R,G,B 

under the simplifying condition that the pairing gap is real, i.e. A* = A. It is then 
straightforward to write down the Heisenberg equations of motion of the field operators: 

d ~ ( h 2 _., . _ L 

(5) 
(6) 
(7) 



in—ipG 
at 



h 2 

2m ' 
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which are coupled by the presence of the same chemical potential /i in the three equations. 
We now use the Bogoliubov-Valatin representation of the field operator ip a (r,t) in terms of 
the anticommuting quasi-particle Bogoliubov operators b ka : 

Mr, t) = J2 (u k e^ r -^ b kR - v k e-^ r -^ b+ G ) , (8) 

k 

Mr,t) = ^(w fc e* (k - r - Wfet) S kG + ^e- i ( k - r -^*)S+ R ) , (9) 

k 

$ B (T,t) = J^^^^B , (10) 
k 

where u k and v k are such that u\ + v\ = 1. After inserting these expressions into the 
Heisenberg equations of motion of the field operators we get 

£ fc = Ml k = — — - // , (11) 
2m 



E k = frw k = Ve| + A 2 , (12) 

and also 



By imposing the following ground-state averages 

(K B b k , B ) = 0(-&) 5 k ,k' , = (bt G k G ) = ®{~E k ) 5 k ,k' , (14) 

with Q(x) the Heaviside step function, the number equation (T5|) gives 

N = N R + N G + N B , (15) 

where 

N R = N G = l -Y< v l ( 16 ) 

k 

and 



k 

Similarly, the gap equation (J3]) gives 



2m 



5 V ^ 2E k 

The chemical potential /i and the gap energy A are obtained by solving equations (Tl5|) and 
f fTSj) . In the continuum limit, due to the choice of a contact potential, the gap equation ( H~8|) 



diverges in the ultraviolet. This divergence is linear in three dimensions and logarithmic in 
two dimensions. We shall face this problem in the next two sections. 

Another interesting quantity is the the number of red-green pairs in the lowest state, i.e. 



the condensate number of red-green pairs, that is given by 
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30. 



32 



33] 



iV n 



d 3 r, d 3 r 2 |<^ G (r0^(r 2 )}| 2 = 5>l 



(19) 



In the last years two experimental groups 



26M28I] have analyzed the condensate fraction 
of three-dimansional ultra-cold two-hyperfine-component Fermi vapors of 6 Li atoms in the 
crossover from the Bardeen-Cooper-Schrieffer (BCS) state of Cooper Fermi pairs to the Bose- 
Einstein condensate (BEC) of molecular dimers. These experiments are in quite good agree- 



ment with mean-field theoretical predictions 



29 



30| and Monte-Carlo simulations |3ll 



temperature, while at finite temperature beyond-mean-field corrections are needed 



at zero 



33]. 



Here we show how to calculate the condensate fraction Nq/N for the three-component Fermi 



gas at zero temperature 



251 ] in three 29J and two |34( dimensions. Finally, we calculate the 



first sound velocity c, of the three-component system by using the zero-temperature ther- 



modynamic formula 



35] 



(20) 



n dfx 
m dn 

where /i is the chemical potential of the Fermi gas and n the total density. The sound 
velocity c s , which is the Nambu-Goldstone mode of pairing breaking of SU(3) symmetry, 
has been previously analyzed by He, Jin and Zhuang 18j] in the three dimensional case. 
Here we study c s in the two dimensional case too. 



III. THREE DIMENSIONAL CASE 



In three dimensions a suitable regularization |2l|, [23j of the gap equation (fTS]) is obtained 
by introducing the inter- atomic scattering length ap via the equation 

1 \ ^ m 



1 ml x ^ m 

g Anh 2 a F V ^ h 2 k 2 



and then subtracting this equation from the gap equation ([18]) . In this way one obtains the 
three-dimensional regularized gap equation 

m 



4nh 2 a F V \2E k h 2 k 2 



1 x ^ / 1 m \ 



In the three-dimensional continuum limit — > V/ {2ir) 3 j <i 3 k — > V/ (27r 2 ) J k 2 dk from the 
number equation f|T5|) with f|T6|) and f|T7|) we find the total number density as 



iV . , 

n = — = n R + n G + n B , (23) 



with 

1 (2m) 3 / 2 

^ = " G = 2^r 

and 



A 3 / 2 / 2 (£) , (24) 



3 2vr 2 ft 3 

The renormalized gap equation (1221) becomes instead 



«^4^lC^f 3/2 eM. (25) 



1 2(2^/^,^ 



op nh \A 

where fc F = (67riV/(3lO) 1/3 = (2tt 2 n) 1 / 3 is the Fermi wave number. Here h(x) and h(x) are 
the two monotonic functions 

vV - xf + 1 I/ 2 



= / [ 7=1= ? ~ -, )dy, (27) 



o 



'• ( *>=' s2 1 1 - TW^Wtt I* • (28) 







which can be expressed in terms of elliptic integrals, as shown by Marini, Pistolesi and 
Strinati [23] . In a similar way we get the condensate density of the red-green pair as 



This equation and the gap equation (1261) are the same of the two-component superfluid fermi 
gas (see [29]) but the number equation (|T5l) . with ( IT6]) and (jTTI) . is clearly different. Note 
that all the relevant quantities can be expressed in terms of the ratio 

x = £ , (30) 

where x e] — oo,oo[. In this way the scaled energy gap A/ep and the scaled chemical 
potential /V e F rea d 

A X ° (31) 



e F i 2 {x ) + lx 3 /2 e{x ) 
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- = TTh > ( 32 ) 

where e F = h 2 k 2 F /(2m) = (2vr 2 n) 2/3 ft 2 /(2m) is the Fermi energy of the 3D ideal three- 
component Fermi gas with total density n. The fraction of red fermions, which is equal to 
the fraction of green fermions, is given by 

UR = nG = h(xo) , . 

n n 2I 2 (x ) + fzf 0(x o ) ' 

while the fraction of blue fermions reads 

M =1 M£|> . (34 ) 

The fraction of condensed red-green pairs is instead 



n 7i \Jx Q + ^/l + xl 

(35) 



n 8V2/ 2 (x ) + |xf 9(x ) ' 
Finally, the adimensional interaction strength of the BCS-BEC crossover is given by 

= J_ = _2 IiM 

k F a F 7\ ( ' , \ , \ -ill r\i \\^ 3 



(j 2 (xo) + l^o 2 ©(^o) 



(36) 



We can use these parametric formulas of xq to plot the density fractions as a function of the 
scaled interaction strength y. 

In the upper panel of Fig. [1] we plot the energy gap A (in units of the Fermi energy 
e F ) as a function of scaled interaction strength y = l/^kpCLF)- As expected the gap A is 
exponentially small in the BCS region (y <C —1), it becomes of the order of the Fermi energy 
ep at unitarity (y = 0), and then it inceases in the BEC region (t/ > 1). In the lower panel of 
Fig. Q]we show instead the scaled chemical potential /i/e^ as a function of scaled interaction 
strength y = l/^kpap)- In the BCS region (y —1) the chemical potential fi is positive 
and practically equal to the Fermi energy e F of the ideal gas; at unitarity (y = 0) the /i is 
still positive but close to zero; it becomes equal to zero at y ~ 0.6 and then diminishes as 
— y 2 (half the binding energy of the formed dimers). 

In the upper panel of Fig. |5] we plot the fraction of red fermions n^jn (solid line) 
and the fraction of blue fermions ub/u (dashed line) as a function of scaled interaction 
strength y = l^kpap). The behavior of nc/n is not shown because it is exactly the 
same of n R /n. The figure shows that in the deep BCS regime (y <^ —1) the system has 



y 



FIG. 1: Ultracold fermions in three-dimensions. Upper panel: scaled energy gap A/ep as a function 
of scaled interaction strength y = l/(kpap). Lower panel: scaled chemical potential n/ep as a 
function of scaled interaction strength y = l/(kpap). Note that kp = (27r 2 n) 1//3 is the Fermi 
wave number and ep = h 2 kp/(2m) = (2ir 2 n) 2j/3 h 2 /(2m) is the Fermi energy of the 3D ideal three- 
component Fermi gas with total 3D density n. 

nn/n = nc/n = ns/n = 1/3. By increasing y the fraction of red and green fermions 
increases while the fraction of blue fermions decreases. At y ~ 0.6, where fi = 0, the fraction 
of blue fermions becomes zero, i.e. ns/n = and consequently ur/u = ric/n = 1/2. For 
larger values of y there are only the paired red and green particles. This behavior is fully 



consistent with the findings of Ozawa and Baym [19]. In the lower panel of FigJT] it is 
shown the plot of the condensate fraction n /(n/2) of red-green pairs through the BCS- 
BEC crossover as a function of the Fermi-gas parameter y = l/{kpap). The figure shows 
that a large condensate fraction builds up in the BCS side already before the unitarity limit 
(y = 0), and that on the BEC side (y 3> 1) it rapidly converges to one. 

As previously stressed, by using Eq. fl20l one can obtain the first sound velocity. In 
particular, we have found that ji = ep F(y), where F(y) is the numerical function plotted 
in the lower panel of Fig. [TJ It is then straightforward to show that 



c s = ^F{y)- l -yF'(y). (37) 
By using this formula we plot in Fig. |3]the scaled sound velocity c s /vp as a function of scaled 
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FIG. 2: Ultracold fermions in three-dimensions. Upper panel: fraction of red fermions nn/n (solid 
line) and fraction of blue fermions ub/u (dashed line) as a function of scaled interaction strength 
y = l/ikpCLF)- Lower panel: condensed fraction of red-green particles uq/ti as a function of scaled 
interaction strength y = l/{kpa,F). Units as in Fig. 1. 

interaction strength y = I/^A^of). The curve shows that c s jvp decreases by increasing y 
and it shows a knee at y c± 0.6, where the chemical potential changes sign. 



IV. TWO DIMENSIONAL CASE 



A two-dimensional Fermi gas can be obtained by imposing a very strong confinement 
along one of the three spatial directions. In practice, the potential energy E p of this strong 
external confinement must be much larger than the total chemical potential fi^p, of the 
fermionic system: fi^p, <C 2Ep [37]. Contrary to the three-dimensional case, in two dimen- 



sions quite genera 
g between atoms 



ly a bound-state energy e# exists for any value of the interaction strength 



22 



231 ] . For the contact potential the bound-state equation is 



1 = 1 E 



k 2m ~ B 



(38) 
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FIG. 3: Ultracold fermions in three-dimensions. Scaled sound velocity c s jvp as a function of scaled 
interaction strength y = l/ikpap). Units as in Fig. 1, with vf = frkp/m the Fermi velocity of the 
3D ideal three-component Fermi gas. 



and then subtracting this equation from 



dimensional regularized gap equation 22 



he gap equation (1T8]) one obtains the two- 



t + eR 2E k 

k V 2m ~ CB K 







(39) 



Note that, for a 2D inter- atomic potential described by a 2D circularly symmet- 
ric well of radius R and depth U , the bound-state energy e B is given by ee — 
h 2 /(2mR 2 ) exp (-2H 2 / (mU R 2 )) with E/ i2g -> |36|. 

In the two-dimensional continuum limit — > V/(2-k) 2 J d 2 \t — > V/(2tt) J kdk, the Eq. 
039]) gives 

(40) 



e B = A 



1 + 



A 2 A 

Note that here V is the 2D volume of the gas, i.e. an area. Instead, the number equation 
(|T5|) with ( JT6]) and ( TTT|) gives the total number density as 

N 



n = — = n R + n G + n B 



(41) 



where V is a two-dimensional volume (i.e. an area), the red and green densities are 



^ = ^ = i^ A U + V 1 + A^ 



(42) 
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while the blue density reads 



■2nh 2 

Finally, the condensate density of red-green pairs is given by 

1 / m 



vn 



(^)A(| + arctan(£)). (44) 



Also in this two-dimensional case all the relevant quantities can be expressed in terms of 
the ratio x$ = fi/A, where xq e] — oo, oo[. In particular, the scaled pairing gap is given by 

A 3 



e F Xq + a/1 + Xq + X Q(x ) 

while the scaled chemical potential reads 

fi 3x 



(45) 



(46) 



£F Xq + a/1 + X 2 Q + X Q(x ) 

where the two-dimensional Fermi energy = h 2 k F /(2m) of the 2D ideal three-component 
Fermi gas with 2D total density n is given by ep = 7rh 2 n/m with kp = (47TO/3) 1 ' 2 the Fermi 
wave number. 

The fraction of red fermions, which is equal to the fraction of green fermions, is given by 

Ur = UG = Xq + a/1 + X% 

n n 2[x + ^l + x 2 + x Q{x )} ' 



the fraction of blue fermions is 



^ = 1-2^, (48) 
n n 



and the condensate fraction is 

n f + arctan (x ) 



(49) 



n 4[x + a/1 + 4 + x 9(x )] 
It is convenient to express the bound-state energy e# in terms of the Fermi energy ep. In 
this way we find 

€ -» = 3 V^o-^ > (50 ) 

e F X + a/1 + x\ + X 9{x ) 

We can now use these parametric formulas of xq to plot the fractions as a function of the 
scaled bound-state energy e^/e^- 

In the upper panel of Fig. H]we plot the scaled energy gap A/ep as a function of scaled 
binding energy es/ep. The gap A is extremely small in the "BCS region" (e#/ep <C 1/2), 
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£ B /£ F 

FIG. 4: Ultracold fermions in two-dimensions. Upper panel: scaled energy gap A/cf as a function 
of scaled bound-state energy e^/e^. Lower panel: scaled chemical potential (J,/ep as a function 
of scaled bound-state energy €b/cf- Note that kp = (4-7rn/3) 1//2 is the Fermi wave number and 
ep = h 2 k F /(2m) = (4irn/3)h 2 /(2m) is the Fermi energy of the 2D ideal three-component Fermi 
gas with total 2D density n. 

it becomes of the order of the Fermi energy ep at tsjtF = 1/2, and then it inceases in the 
"BEC region" (e^/e^ 3> 1/2). In the lower panel of Fig. 0] we show instead the scaled 
chemical potential fi/ep as a function of scaled binding energy e^/e^. In the BCS region 
(cb/cf <C 1/2) the chemical potential /i is positive and decreases as /i = tF — 6b/3; /i becomes 
equal to zero at €b/cf = 3 and then it further decreases linearly as [i = (3€f — es)/2. 

In the upper panel of Fig. Owe plot the fraction of red fermions n#/n (solid line) and the 
fraction of blue fermions ns/n (dashed line) as a function of scaled bound-state energy es/^F- 
The behavior of nc/n is not shown because it is exactly the same of n^/n. The figure shows 
that in the deep BCS regime (e^/e^ <S 1) the system has n R /n = n G /n = N B /n = 1/3. By 
increasing €b/^f the fraction of red and green fermions increases while the fraction of blue 
fermions decreases. At ts/^F = 3, where /x = 0, the fraction of blue fermions becomes zero. 
For larger values of €b/€f there are only the paired red and green particles. This behavior 
is quite similar to the one of the three-dimensional case; the main difference is due to the 
fact that here the curves are linear. In the lower panel of FigfS] it is shown the condensate 
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FIG. 5: Ultracold fermions in two-dimensions. Upper panel: fraction of red fermions n/j/n (solid 
line) and fraction of blue fermions rig/n (dashed line) as a function of scaled bound-state energy 
€b/cf- Lower panel: condensed fraction of red-green particles no/n as a function of scaled bound- 
state energy eb/^f- Units as in Fig. 4. 

fraction n /(n/2) of red-green pairs. In the weakly-bound BCS regime (e^/e^ <C 1) the 
condensed fraction no/n goes to zero, while in the strongly-bound BEC regime (e^/e^ 3> 1) 
the condensed fraction no/n goes to 1/2, i.e. all the red-green Fermi pairs belong to the 
Bose-Einstein condensate. Notice that the condensate fraction is zero when the bound-state 
energy 6b is zero. For small values of e^/e^ the condensed fraction has a very fast grow but 
then it reaches the asymptotic value 1/2 very slowly. 

Also in 2D, by using Eq. (120]) one can obtain the first sound velocity. We have found 
that fi = €f G(€b/€f), where G(€b/€f) is the numerical function plotted in the lower panel 
of Fig. HI It is then straightforward to show that 



v^V e F e F K e F ' 



(51) 



By using this formula we plot in Fig. [6] the scaled sound velocity c s /vf as a function of 
scaled binding energy e^/e^. The curve shows that c s /vf is constant, i.e. c s /vp = l/v^2, 
by increasing e^/e^ up to ee/e^ = 3 where the chemical potential becomes equal to zero. 
For a larger value of e B /ep the sound velocity c s jumps to a larger constant value, i.e. 
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FIG. 6: Ultracold fermions in two-dimensions. Scaled sound velocity c s /vf as a function of scaled 
bound-state energy e_e/e_p. Units as in Fig. 4, with vf = hkp/m the Fermi velocity of the 2D ideal 
three-component Fermi gas. 

c s / v f — V / 3/2. This kink in the first sound velocity c., is reminescent of the jumps seen with 



repulsive fermions with reduced dimensionalities 



23, 



and with dipolar interaction [39] . 



V. CONCLUSIONS 



We have investigated a uniform three-component ultracold fermions by increasing the 
SU(3) invariant attractive interaction. We have considered the symmetry breaking of the 
SU(3) symmetry due to the formation of Cooper pairs both in the three-dimensional case 
and in the two-dimensional one. We have obtained explicit formulas and plots for energy 
gap, chemical potential, number densities, condensate density, population imbalance and 
first sound velocity in the full BCS-BEC crossover. In our calculations we have used the 
zero-temperature mean-field extended BCS theory, which is expected to give reliable results 
apart in the deep BEC regime 3lM33|. Our results are of interest for next future experiments 
with degenerate gases made of alkali-metal or alkaline-earth atoms. As stressed in the 
introduction, SU(N) invariant interactions can be experimentaly obtained by using these 
atomic species [9|, [lOj, ll3|. The problem of unequal couplings, and also that of a fixed 



number of atoms for each component, is clearly of big interest too, and its analysis can be 
afforded by including more than one order parameter 40]. 
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There are other interesting open problems about superfluid ultracold atoms we want 
to face in the next future. In particular, we plan to investigate quasi one-dimensional 
and quasi two-dimensional Bose-Einstein condensates in nonlinear lattices (i.e. with space- 
dependent interaction st reng th) 41] . Moreover, we want to analyze the signatures of classical 



and quantum chaos 42h46| with Bose-Einstein condensates in single-well and double-well 



47 



491 ] . Finally, we aim to calculate 



configurations, and also in the presence of vortices 
analytically the coupling tunneling energy of bosons by means of the WKB semiclassical 
quantization 



50 



53[ and comparing it with the numerical results of the Gross-Pitaevskii 



equation. 
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